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The grand partition function of a model of confined quarks is exactly calculated at arbitrary 
temperatures and quark chemical potentials. The model is inspired by a softly BRST-broken version 
of QCD and possesses a quark mass function compatible with nonperturbative analyses of lattice 
simulations and Dyson-Schwinger equations. Even though the model is defined at tree level, we 
show that it produces a nontrivial and stable thermodynamic behaviour at any temperature or 
chemical potential. Results for the pressure, the entropy and the trace anomaly as a function of the 
temperature are qualitatively compatible with the effect of nonperturbative interactions as observed 
in lattice simulations. The finite density thermodynamics is also shown to contain nontrivial features, 
being far away from an ideal gas picture. 


I. INTRODUCTION 

The study of the infrared behavior of non-Abelian gauge theories, like Quantum Chromodynamics, presents a 
remarkable task for theoretical physics. The strongly-coupled nonperturbative nature of infrared phenomena greatly 
undermines our ability to access the physical properties of the theory in an analytical way. One general question 
raised in this context concerns the description of the degrees of freedom in a confining theory. Confinement of the 
perturbative excitations associated with the fundamental fields means, among other features, that these degrees of 
freedom do not appear in the physical spectrum of the theory as asymptotic states. The physical spectrum in this 
case is given by bound states associated with composite fields of the theory. The fundamental problem is thus to 
understand how the perturbative degrees of freedom disappear from the spectrum and how this is reflected in an 
action formulation written in terms of them. 

A very important tool in identifying the physical space of states in a gauge theory is the BRST symmetry [312 ■ In 
a BRST invariant theory, the physical states are defined by the cohomology of the nilpotent BRST operator that acts 
on the fundamental fields. However, this is strictly valid only at the perturbative level. In fact, the whole construction 
of the BRST formalism in a gauge theory can be traced back to the Faddeev-Popov procedure of gauge fixing. Since it 
has long been known that this procedure ceases to be valid at the nonperturbative level [3] , there is no strict necessity 
to expect that the BRST symmetry survives in the infrared regime. Indeed, there is strong evidence coming from 
lattice studies that the BRST symmetry is broken at low energies |3]- This is in tune with the expectations for a 
confining theory: the survival of the original BRST symmetry would allow for a definition of the physical spectrum, 
by its cohomology, in terms of the fundamental fields, thus contradicting confinement. 

It is therefore expected that the description of the original perturbative degrees of freedom is modified in the 
infrared. This seems to be exactly what happens in pure gauge theories according to lattice results concerning the 
Yang-Mills gauge field and ghost propagators in the Landau gauge m\- positivity violation is observed for the gluon 
two-point function. It is well-known 011] that the two-point function of a physical field is directly related to the 
probability of propagation of the associated asymptotic state. Therefore, positivity violation, as observed for lattice 
gluons, hinders the probabilistic interpretation of the two-point function in terms of propagating particles, being 
consistent with confinement in the sense of absence of colored states from the spectrum. This is the confinement 
criterion we adopt (more discussion can be found in e.g. the reviews [Mill]). 

The numerical data also point to gluon propagator that tends to a constant value in the infrared, while the ghost 
propagator has the behavior of a free particle. In the context of the Schwinger-Dyson equations, this is called the 
decoupling solution m- The theoretical description of these results in the continuum has been proposed to be 
provided by an effective action in the so-called refined Gribov-Zwanziger (RGZ) formalism [T3Ill5j . The RGZ action 
is renormalizable and displays a softly broken BRST symmetry. One very important feature of this formalism is the 
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guarantee that the quantum ultraviolet behavior of the theory is under control and the - very successful - original 
perturbative Yang-Mills framework is recovered at high energies. 

The RGZ action describes a pure gluonic theory and its construction has the geometrical guidance of the original 
Gribov observations about the configuration space of the gauge fields. Its action effectively describes the restriction 
of the path integral measure to a region without the infinitesimal Gribov gauge copies, while also taking into account 
the formation of condensates. The resulting action has tree-level propagators that describe very well the lattice data 
at zero and finite temperature naEzi and can be used to obtain consistent estimates of the masses of glueball states 
with different quantum numbers [MIS]. 

An interesting development is that lattice data on propagators suggest a similar behavior for the gauge-interacting 
matter fields, displaying positivity violation [20], as observed for the gluon fields. Furthermore, the theoretical 
description of these degrees of freedom in the infrared has been conjectured to have the same structure as for the 
gluon, with an action whose tree-level propagators reproduce very well the lattice data, also displaying softly broken 
BRST symmetry [^. It should be noticed, however, that in the case of matter fields a geometrical interpretation of 
the modified action as a restriction in the configuration space of the fields is lacking. 

This set of affairs led us to expect that a possible signature of confinement would be the soft breaking of the 
BRST symmetry. In fact this relation has been explored in many works |22H2H]- The idea is that confined degrees 
of freedom are described by an effective action with a BRST-breaking soft term. Interestingly, the breaking of the 
BRST invariance may have a connection with chiral symmetry breaking, since the currently known formulations of 
soft BRST breaking in the quark sector always imply chiral symmetry breaking |3D| and thus dynamical mass 
generation in the deep infrared seems to be intimately related to BRST breaking m- In this sense, one may hope 
that the investigation of the nonperturbative breaking of the BRST symmetry may shed some light on a long-standing 
puzzle in the QCD phase structure: the apparent and unexplained link between confinement and chiral symmetry 
breaking. 

Several effective models of QGD present in the literature address the relation between chiral symmetry breaking 
and confinement. Two important prototype models extensively explored to address the question of chiral symmetry 
breaking in the strong interactions are the Nambu-Jona-Lasinio (NJL) and the Quark Meson (QM) models |5T1132j . 
generically called chiral models. On one hand, chiral symmetry breaking and restoration is, in chiral models, a result 
of the interaction of the quark sector either to itself (in a four-fermion vertex, as in the NJL model) or to meson 
fields (as in the QM model). On the other hand, as an attempt to take into account the physics of confinement, 
one can introduce a Polyakov loop field coupled to quarks [551 - 155] . Within the Polyakov loop extended NJL (PNJL) 
and QM (PQM) models, the Polyakov loop field effectively encodes the gauge field as an external background that 
interacts with the quark degrees of freedom, so that the deconfinement transition may be studied in this context as 
a (static) Landau-Ginzburg system coupled to the chiral model. In spite of their interaction with the Polyakov loop, 
the quark degrees of freedom in both the PNJL and PQM models physically correspond to asymptotic states, not 
to confined ones. This is so because their two-point functions only have poles on the real axis and thus respect the 
Osterwalder-Schrader’s axiom of reflection positivity [39], guaranteeing the establishment of a standard probabilistic 
description of propagation of asymptotic quarks, in contradiction with real-world QCD. 

A possible way to extend the NJL model is to include a nonlocal current-current interaction kernel [5DIH5] . As a 
result, complex quark masses may appear, indicating confinement in the sense we discussed before. However, due to 
its four-fermion interaction, the nonlocal NJL model is not renormalizable, similarly to its local version. Besides, it 
also leads to unstable thermodynamic behavior [42ll44j . 

In the last years, some works have explored the thermodynamic properties of softly BRST broken pure gauge 
systems (see, for example, |45H50j L In this paper, we explore thermodynamic properties of a model of quarks with 
soft BRST breaking. As discussed above, this model is expected to describe the infrared properties of confined 
(positivity-violating) quarks, while keeping compatibility with ultraviolet QCD properties. Indeed, the analytical 
propagator of the model fits well the available lattice data [50] and the model has been proven to be renormalizable 
[29| . reducing to perturbative quarks in the ultraviolet regime. Our goal here is to show not only that soft BRST 
breaking in the quark sector implies a well-defined macroscopic behaviour, but also that the tree level model is 
capable of predicting nontrivial features, being in general qualitatively compatible with the effect of nonperturbative 
interactions as observed in lattice data. Furthermore, the instabilities present in the nonlocal NJL model are absent 
from our setup. 

This work is organized as follows. Section II presents the action dehning the model. In Section III the exact 
partition function is computed within the imaginary-time formalism. In section IV we obtain the thermodynamic 
quantities and present our results for the distribution function, the pressure, the entropy, and the trace anomaly. Our 
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conclusions are discussed in Section V. 


II. THE ACTION OF THE QUARK MODEL WITH SOFT BRST BREAKING 

Following the discussion of |29j , let us briefly review how a model with soft breaking of the BRST symmetry can be 
obtained from an extension of the QCD lagrangian. We start from the gauge-fixed lagrangian density in an euclidean 
space of dimension 4, 


^QCd = + c + ed^Dfc\ (1) 

where — d^A^^ + gA^^A^ is the field strength and D'-^ = 5*-^9^ — ig{T‘^y^A’^^ is the covariant derivative. 

The indices {a,b,... ,h) = 1,..., — 1 correspond to adjoint representation and = 1,..., W correspond 

to fields in the fundamental representation of SU{Nc). The greek indices (a, /3,...) = 1, 2, 3,4 denote spinor indices 
while ...) = 1,2,3,4 regard euclidean space indices. For more details on notational conventions, the reader is 
referred to Ref. [5T] . 

One very important feature of Q is its invariance under BRST transformations mm, 

sA- = 

syi = -igc%T<^ryi 
= -ig^ic^iT^r 

sc“ = Igr^^c^c- 

= ib‘^ 

sb’^ = 0 ( 2 ) 

which is crucial for its multiplicative renormalizability. As an intermediate step in the definition of the model, let 
us introduce two BRST doublets (^*,0*) and ( 77 *, A®) of auxiliary helds (and their hermitian conjugate fields) that 
transform as 


The addition of the BRST invariant action 


= C, sOl = 0 
SVa = K, sA), = 0 . 


S'jA = s 


J d^x - Cvl,)] 


[-AL9^CL - (ALCL + CAj. + 


(3) 


(4) 


does not change the physical content of the theory. This can be easily checked by integrating the auxiliary fields, 
which trivially gives 1. However, if one also adds to the action the coupling 


Sm = J d^x + ca) - M^CKr^ + CAD] 


(5) 


between the auxiliary fields and the matter fields, the resulting action Stot = Sqcd + -I- Sm is no longer BRST 

invariant. However, as shown in [29j . this does not destroy the renormalizability of the theory once it corresponds to 
a soft breaking of the BRST symmetry. In other words, the theory modified by Sm is equivalent to actual QCD in 
the ultraviolet. 

Of course, the addition of ([^ leads to changes of the theory in the infrared with respect to actual QCD and one 
might ask to what extent such a theory can correctly describe the infrared physics of QCD. Interestingly enough, 
BRST soft breaking terms can be dynamically generated. For example, gluon condensation may be incorporated in 
the action and serve as a starting point for an effective theory containing nonperturbative physics |22j . Indeed, the 
presence of such condensates leads to a nontrivial infrared behavior of gluon and ghost propagators, which can be also 
observed in lattice simulations mm- In the quark sector of the theory, the BRST breaking term (§ is conjectured to 
be a consequence of the condensation of the qq operator. 
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The physical meaning of the action Stot becomes clearer after integration in the auxiliary fields 
and (A, A). As a result, due to the BRST breaking term Sm, one finds a nonlocal action 


Snl = / <rx V’a 




2MfM2 

—9^ + m? 


+ Too 




( 6 ) 


for the quark sector, plus the minimal coupling to the gauge sector. As a first approximation in our model, we neglect 
the interaction between the gauge and the matter sectors, except through the nonperturbative term Sm- Notice that 
the usual local free fermion case can be recovered for 2 M 1 M 2 = M 3 —)■ 0 . 

In the vacuum, the tree level fermion propagator reads (with omission of color and Dirac indices) 


S{P) 


+ Mpjp) 

+ M§{p) 


where we defined the vacuum mass function 


Mo(p) 


2 MIM 2 

+ vn? 


+ mo- 


(7) 

( 8 ) 


Such a quark mass function is compatible with lattice QCD calculations at zero temperature [20] as well as with an 
analysis of the Dyson-Schwinger equations of the quark propagator |SS]. The parameters M 3 = 2MfM2, rn^ and mg 
of the mass function ^ can be well fitted to the lattice results of [20], giving M 3 = 0.196GeV^, = 0.639GeV^, and 

Too = 0.014GeV with x^/dof = 1.18 [30]. Notice also that the propagator Q has two complex conjugate poles and 
a real pole. This can be interpreted as a sign of confinement, a property expected from strongly interacting particles 
at zero temperature. 

Once the action of the quark sector 


Sq= So,pmo]tpp + S^x + Sm- 


(9) 


is quadratic in the fermion fields, it is possible to calculate the partition function exactly for arbitrary temperatures 
and chemical potentials. In the next section, we will calculate the grand canonical partition function of the theory 
defined by 


III. THE PARTITION FUNCTION : AN EXACT COMPUTATION 

In order to calculate the partition function of the theory ® , we shall use standard techniques of hnite temperature 
field theory in the Matsubara imaginary time formalism [51ll56j . At finite temperature, the 0-direction is compactified 
(M —^ = [0,/3)), with P = 1/T and the fermion helds obey anti-periodic boundary conditions in imaginary time. 

As a consequence, the fermion fields can be written in terms of their Fourier transforms as 

^/>(x4,x)=^y (10) 

and 

v;(x4,x)=^ y (11) 

where a;„ = {2n + \)'kT are Matsubara frequencies (n G Z). 

The presence of a net quark charge requires a nonzero quark chemical potential p. Indeed, a Noether charge can 
be associated with a global 17(1) transformation of the type 

pj -G ( 12 ) 

with a S K. It is crucial to notice that the action © is invariant under this transformation only if all the auxiliary 
fields also carry this same charge. The corresponding Noether current is then used to impose the grand canonical 
constraint on the hamiltonian of the theory. The calculation is straightforward but it must be carefully performed. 
The details of the derivation of the in-medium effective action for the quark helds are given in the Appendix 
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After imposing the grand canonical constraint and integrating out the auxiliary fields and their respective momenta, 
the grand canonical partition function can then be written in the nonlocal form 


Z{T,^) = Tr exp —j3(^H — ^N^ 


(13) 


= / [Dijj] [Dijj] exp 




where, omitting color and spinor indices. 


= Ip 


74 (i94 - m) - *7 • V + 


M 3 

—{84 — ~ 


+ Too 


Y- 


(14) 


That is, our calculation explicitly showed that the inclusion of the chemical potential simply amounts to the shift 
84 —> 84 — fi. It follows from (101 and 0 that 


^ni - 


(Yp 

(^ 


V’(-w„, -p) [P (-74(iw„ + /r) - 7 • p + Mn^pip))] YiuJn, p) 


(15) 


where 


Mn,p{p) ■■= 


M 3 


-{icOn + pY + p2 + 


■ Too- 


(16) 


corresponds to the mass function Q at finite temperature and chemical potential. 

Being quadratic in the fermionic operators, the partition function can be formally integrated immediately, giving 


Z{T,p)= det [-/ 374 (zw„ + ^) - ^7 • p+ /3M„_p(^)] 

P,DJ,C 


(17) 


where the (full) determinant has to be taken with respect to the momentum (p), Dirac (D), flavor (f), and color (c) 
subspaces. 

Let us follow the same steps as in, e.g., [5S], in order to calculate the grand partition function 0 . The determinant 
in flavor and color spaces are very simple, once all Nf flavors and Nc quark colors are on equal footing. This corresponds 
to taking the determinant in Dirac and momentum subspaces to the power N^Nf. Furthermore, the Dirac determinant 
can be straightforwardly computed to be 

det [-/ 3 j 4 iiuJn + p)- Pj-p + (3Mn^p{p)] = /3'^ [p^ + M^p{p) - (jw„ + pY]'^ ■ (18) 

Before proceeding to the calculation of the determinant in momentum space, it is convenient to consider the 
logarithm of the partition function and use the operator identity log det A = Tr log A where, in momentum space, 
A = [p^ + M^pip) + (w„ — ipY] ■ The trace in momentum space correponds to summing over all Matsubara 

frequencies and all momenta. In the thermodynamic limit (V —>■ 00 ), 

P 


As a result, we find 


log Z{T, p) = 2VN,Nf E / ^ log [p' + + pf] 


( 20 ) 


Notice that the sum (20) is real although its integrand is not. Using that ujq = —W-i, wi = —uj -2 etc., it is not 


difficult to write (20) explicitly as a sum of real numbers by splitting the sum of Matsubara frequencies for n > 0 and 


n < — 1. In order to continue the calculation, it is best to split Eq. (20) in two factors, adding and subtracting the 
p = 0 contribution, so that 


log Z{T, p) 

2pVNcNf 


^ log [p 2 + Mlp{ 0 ) + } + ^ log ■ 


+ MYp{p) - jiuin + pY 

p2 + MYp(O) + Lo^ 
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_ iogz(r,o) iogz(^)(r,Ai) 

2/3VNcNf 2pVNcNf '■ 

where we used the standard sum-integral notation, 


« oo 


tfp 

{ 2 tt)^ 


(•••)• 


( 21 ) 


( 22 ) 


One interesting feature of the decomposition (21) is that only the term log ^(T, 0) contains any divergencies. They 


are the same as those of the local theory at T = 0 and therefore can be straightforwardly subtracted. The /r 7 ^ 0 
contribution is finite, as expected from the argument that the introduction of a chemical potential should not bring 
any new divergencies to the theory. For this reason, it can be calculated numerically and also has a closed expression 
in terms of elementary integrals in the limit of zero temperature. 

For clarity, let us now separately discuss the calculations of the /r = 0 and p ^ 0 contributions. 


A. The p. = 0 contribution 


In the standard calculation of the partition function of local free quarks of mass toq at ^ = 0 (see, e.g., [5T1 IbH] !. 
one is lead to the calculation of the sum-integral 


log ZiociT, 0 ) = 2 N,NfPV ^ log [luI + } , 


(23) 


where uj = + TOq. After subtractions of T- and ^-independent (infinite) constants, one arrives at the hnal 

expression 

d?p 


logZioc{T,0)=2N,NfPV J ■^[uj + 2T\og{l + e-^‘^)]. 

It is quite straightforward to show that the calculation of 

log ZiT, 0 ) = 2l3VN,Nf ^ log {/3^ [p^ + ^( 0 ) + u;^] } 


(24) 


(25) 

can be reduced to the sum of four terms which are each formally identical to ( |23[ ). Indeed, using Eq. ( [I^ , it is 
possible to write the argument of the logarithm in (|25|), Arg, as a ratio of two polynomials, 


^ [(P^ + ml+ w2)(p2 + m^ + a;2)2 -p 2M3mo{p^ + + ujD + M|] _ 

■^'^9 •— , ..9\9 = P 




(p 2 -I- -I- a; 2)2 


(26) 


(p 2 + m'^ + a; 2)2 

where the third degree polynomial PaioJ^) can be factored out as 

= (ujI + pDicl + + pI), (27) 

with —ipf (i = 1,2,3) the three roots of P3. In general, the roots of a third degree polynomial of real coefficients 
are one real number and a couple of complex conjugate numbers. This is indeed the case for our parameter set, as 
should be evident from the presence of complex poles in the propagator, i.e., complex masses. The thre e roots can be 
explicitly calculated as functions of p^ and the parameters M 3 , mp, and mg. Substituting back in (251, we have 

log Z{T, 0) = 2N,NfpV ^ IE log K + - 2 log [P^ + p^ + m^)] | . (28) 


As advertised, each term has the same structure as (|23|). Therefore, it follows straightforwardly that 


log Z{T, 0) = log Zo -h ANcNfV 


d^p 


log 


(1 


=-/3vi 


)(1 


=-/3V2 


)(1 


where 


(27r)3 

dPp 


log Zo = 2NcNfl3V j +ip 2 + iP3- 2 <^o) 


(29) 


(30) 


is the pure vacuum contribution and vjq = \/p^ -|- rrp. We choose a normalization such that log Z(0, 0) = 0, i.e., 
logZg = 0. Finally, notice that although two of the (pi are complex, their imaginary parts cancel in (29), so that the 
final result is real, as it had to be. 
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B. The finite fi contribution and its zero temperature limit 


The /i ^ 0 contribution in (211, 


log (T, = 2pVN,Nf ^ log . 






M2 

n,p 


( 0 ) 


■ w" 


(31) 


for the partition function is finite and can be calculated numerically for any finite T and /i. Unfortunately, we could 
not find a closed expression for it, except in the interesting T —> 0 limit, where it equals the full partition function. 
The zero temperature limit is taken most easily by expressing the Matsubara sum in ( |31[ ) as a contour integral in the 
complex plane [56]. It is adequate to choose the rectangular integration path as V =~(—iL + /r + e, +iL + /r + e) U 
{+iL + ^ + e, +iL + p, — e) U {+iL + /r — e, —iL + /i — e) U {—iL + /x — e, —iL + /x + e), with L —)• oo and e —)■ O’*". Defining 


and 


f^p(C) :=P^ + 


/(C) := log ■ 


M, 


n 2 




■ Too 


^^|[(C-/^)^]-(C-/^)M ’ 


we can write the sum over Matsubara frequencies in (|31|) as 


(32) 


(33) 




/(za;„ +f^) = T£ t^^h ^ ^ 




' —200+/i+e 


^/(atanh 


/3(C - m) 


f* —-ioo+zj.—e 


f 200+/i.—e 


i/(a tank 


/3(C - y) 


In the T ^ 0 {/3 ^ oo) limit, 


tanh 


^(C-m) 


^ 1 - 26» [Re(/x - C)] 


so that it follows from (34) that 


^ioo+/j,+e 


T ^ f{iuj„ + ^i)= I ^ /(C) - f 

'ioo+/j.+e ^ Jzc 

^200 + /J, 


—200+^—e 


ZOO+/2—e 




/ 200 + /2 7^ 

;^/(C) 

-200 + /J. 


' —200+/J. 

" de 

27r 


/(t6» + /x) 


(34) 


(35) 


(36) 


The analiticity of / was also used to drop the e term in the integration limits and, in the last step, we made the 
change of variables ^ ^ 9 = —z(C ~ d)- The zero-temperature limit of the complete partition function is thus 

/ 73 /*00 jr\ 

(^ /_ + 

= 2^VN,Nf I ^ g [f[^e + ^) + f{-^e + /x)]. (37) 

It is important to notice that once f{—i9 + fj,) = f*{i6 + fj,), V0 £ M, the partition function ( |3^ is a real quantity, as 
it had to be. 

As a crosscheck, it is possible to show that for the local case, M3 = 0, one simply has p(^) = p2 -|- Wq and thus 


logZz„,(0,/x) =2/3UfV,A^/ 


d^p 

(27r)3 





(38) 


which is the well known result for free local quarks (ST] . 

This ends our calculation of the partition function of the model at finite temperature and quark chemical potential. 
Notice that no approximations were used after the introduction of the model. In the next section, we will present our 
results for thermodynamics quantities directly calculated from the partition function Z (T, p). 


























IV. THERMODYNAMIC QUANTITIES 


Starting from the partition function of the model, 


log ZiT, 

2pVN^Nf 


= 2T 


(Pp 

(27r)3 


log 


(1 


=-/3vi 


)(1 


=-/3V2 


)(1 


o-PV3 


(1 + e-^voy 


log ■ 


+ i^n - ipf 


it is possible to calculate any thermodynamic property of the system. For example, the pressure is given by 

P{T,p) = ^ log Z{T,p), 


(39) 


(40) 


which equals^ minus the thermodynamic potential energy per unit volume, fl{T,p) = —P(T,p). Therefore, other 
quantities may be written as derivatives of the pressure, such as entropy density. 


or the quark number density 


dfl dP 

s{T,p) = - — {T,p) = —{T,p) 




(41) 


(42) 


The energy density can be derived from the other quantities using the thermodynamic identity e = Ts — P + pn, 
while the quark number susceptibility from the second derivative of the pressure. 


dn 


(43) 


In what follows we analyze different regimes of the thermodynamics of the nonlocal quark model. We recall that 
all quantities are calculated from the partition function (39) with Nc = 3 and Nf = 2. 


A. Thermal case 

Let us first address the thermodynamics of a medium of hot confined quarks at zero chemical potential. In this 
regime lattice simulations provide a robust reference for the full QCD case. 

In Fig. &), we show the pressure as a function of the temperature at zero chemical potential. We normalize the 
result by the pressure of free, massless fermions. 


PsB{T,p)=N,Nf 


'7PT^ 

180 
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(44) 


evaluated at /r = 0. For comparison, we also display the pressure of massless quarks subject to a constant negative 
hag pressure psag = —B = —(0.145GeV)'^. 

A smooth but fast rise of the pressure as temperature increases is observed, indicating a drastic change in the 
number of thermal degrees of freedom between low- and high-temperature systems. A similar behavior is seen for 
the thermal crossover in lattice QCD simulations. In contrast to the bag model, our model is capable of describing 
the whole range of temperatures, with no negative pressures attained at low temperatures. This corroborates the 
thermodynamical stability of the model and its consistency. The associated energy and entropy densities at /i = 0, 
normalized by their respective free gas limits, present a similar behavior as functions of temperature, as can be seen 
in Fig. ([^). 

Another thermodynamic observable of interest is the (normalized) trace anomaly or “interaction measure” 

E — 3P 

A{T,p) = ^^. (45) 


^ As long as the normalization condition logZ(0, 0) = 0 holds. 
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FIG. 1: (a) Pressure as a function of temperature for /r = 0, normalized by the free gas limit (44) for nonlocal 
quarks and for the MIT bag model with bag constant B = (O.ldSGeV)"*. (b) Energy and entropy densities, 

normalized by the respective free gas limits. 


It corresponds to the deviation from the tracelessness of a conformal energy-momentum tensor (normalized by T'^), or, 
in other words, it measures how different the nonlocal quark system behaves with respect to an ideal gas of massless 
quarks. The trace anomaly for the nonlocal quarks, plotted in Fig. is clearly different from the bag model result 
and shows a peak around T ~ O.lSGeV. This suggests a smooth transition between two quasi conformal phases at 
low and high temperatures. However, in such a transition, the sound velocity squared. 


c 


2 

S 


dP 


typically presents a minimum around the transition temperature and this is not seen in ([^) . 


(46) 




FIG. 2: (a) Trace anomaly as a function of temperature for zero chemical potential, normalized by T^. A clear 
distinction from the bag model result is seen, (b) Sound velocity squared as a function of temperature. 


Due to asymptotic freedom and chiral symmetry restoration, encoded in the p —>■ oo limit of the mass function (H, it 
is expected that the system approaches the limit of free, massless quarks as temperature increases. This is indeed seen 
from our results for all the thermodynamic quantities we calculated. On the other hand, the non-trivial interactions 
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included in our model via the nonlocal background show up in the results for intermediate temperatures. There is a 
clear change of behavior of all thermodynamic quantities at T ~ O.lSGeV, the typical temperature scale of the chiral or 
deconfined phase crossover. In particular, the trace anomaly shows a peak precisely at this temperature, with a steep 
rise from zero as observed e.g. in lattice QCD simulations and in contrast to QM models. It is interesting that the 
inclusion of a well-chosen Polyakov-loop potential (fitted to thermal lattice data) may provide exactly this steep rise in 
the trace anomaly at low temperatures, which is absent in pure chiral models [38j . In this sense, and for this particular 
observable, our nontrivial background - encoded in the nonlocal quark propagator htted to zero-temperature lattice 
results - seems to play a similar role as the Polyakov-loop potential, further suggesting its connection to confinement 
and nonperturbative physics. 



FIG. 3: High-temperature fit (dotted line) of the trace anomaly as a function of temperature for zero chemical 
potential, normalized by The solid, red line represent our model results. 


In order to further understand the physics contained in our result for temperatures above the peak, we display in 
Figj^a fit of the trace anomaly in the temperature interval T € [300,800] MeV in the following form: 


[A(r, /r — 0)]highT — a -l- ^2 ”1” 


c 

rp4 y 


(47) 


with a = —0.069, b = 0.062 GeV^ and c = 0.0016 GeV^, with the 0{T ^) term dominant over the 0{T ‘^) one in this 
range. 


A fit of the form (471 was also performed for lattice data in Ref. m- There they observe a similar hierarchy of 
contributions, with 6iatt ~ 0.1 GeV^ and ciatt 0.02 GeV^. We note, however, that direct comparison of these results 
should be taken with care, since the system from Ref. m is (2-|-l)-flavour full QCD, while our model describes only 
heavier-than-physical degenerate up and down quarks^, with gluons not explicitly included (their influence appears 
only in the nonperturbative quark dressing). It is nevertheless encouraging that the same structure and hierarchy of 
contributions is found in our confining quark model, in contrast to the bag model or a free massive quark system, as 
will be detailed below. 


The three terms in the fit (47) may be associated with different physical contributions. The constant a is related to 
a logarithmic dependence in the pressure, i.e. ~ logT, that could in principle be mapped to the perturbative result 
(including the running coupling), if the temperature is high enough. The parameter c is originated by a constant 
pressure, which is negative for positive c. This is exactly the type of contribution included in an ad hoc manner in a bag 
model to mimic confinement, and that appears here as a natural consequence of the nontrivial background considered. 
An estimate of the bag constant predicted in the high-temperature region of our model can then be obtained from 
c = AB, yielding B ~ (141 MeV)'^ which nicely falls in the ballpark of the values adopted in the literature of QCD 
models. 

The 0{T~^) term in the trace anomaly, however, is usually absent or very small in a naive bag model. The smallness 
of this term is directly related to the lightness of up and down quarks. Indeed, a massive bag model has a pressure 


^ The light quarks in Ref. [57] also have larger-than-physical masses, corresponding to m-j: = 220 MeV, while the pion in the lattice data 
Eq] for the quark propagator fitted here has m-n- ~ 240 MeV. 
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of the form |56) : 


^bag = - B , (48) 

where m is the current quark mass. The respective trace anomaly is T^A^ag = 4i? + Taking m = 14 MeV, 

the coefficient of is ~ 10“^ GeV^, three orders of magnitude below the value obtained in this ht of our model 
and also in lattice QCD. It becomes clear that the presence of a dominant 0(T“^) contribution may be linked to 
the generation of a large mass scale in the nonperturbative quark sector, as the M3 parameter in our model. It is 
worthwhile noting that even a bag model that considers effective quark masses large enough to give the desired 
contribution would fail to reproduce lattice data for T < 300 MeV, due to the nonperturbative peak structure that is 
developed for lower temperatures. 


B. Cold and dense confined quarks 

Let us now turn our attention to the T —> 0 case and study some thermodynamic quantities as functions of the 
chemical potential /x. This regime of QCD presents a severe Sign Problem which hinders the application of Montecarlo 
simulations, so that model predictions are extremely valuable tools to explore this region of the phase diagram. 



(a) 

FIG. 4: Pressure as a function of the chemical potential for T = 0, normalized by the free gas limit 
PsB{0,f^) = lVcV//x^/127r^. The massless (dotted) and massive (dot-dashed) free gas results are plotted for 

comparison. We use Mthr = 467 MeV. 


In Fig. Q, we show the pressure as a function of the chemical potential, at zero temperature. We compare the 
pressure from the model with that of local massive quarks of mass Mthr = 0.467 GeV, both normalized to the pressure 
of massless quarks. The choice of the mass Mthr is the one that provides a pressure that better compares to our 
model results in order to make explicit the differences and the specific features that are exclusive to the confining 
quark model here studied. 

Since the system is at zero temperature there is no thermal energy to excite particles at low densities, yielding 
a vanishing thermodynamical response. As already observed in the thermal case, the results are fully stable and 
consistent with general physical expectations for a cold and dense system. In particular, the pressure vanishes for 
chemical potentials up to some value fj, ~ Mthr = 0.467 GeV, which is consistent with a dynamically generated scale 
and not directly connected to any specific mass parameter of the model. Starting at this point, the pressure smoothly 
rises until reaching the limiting pressure of local quarks for high chemical potentials. It is not surprising that the 
massless limit is attained at large chemical potentials, given the /i-dependence of the momentum-dependent mass 
function in Eq. (16). 

Besides the pressure, two other relevant observables that can be calculated are the quark density (42), and the quark 
number susceptibility (43), as shown in Figs. ([^) and ([^). The behavior of these observables is reasonably smooth 
and interpolates between that of free massless and massive quarks. The quark susceptibility is a non-monotonic 
function of the chemical potential, having an oscillation inside the window 0.8 GeV ^ /i ^ 1 GeV, resembling the 
result for the massive gas of free quarks, but with a smoothened behavior. 












12 




(a) (b) 


FIG. 5: (a) Quark density as a function of chemical potential /r (with the results normalized by the free massless 
case in the detail) and (b) quark number susceptibility at T —>■ 0. Notice the approach to the free, massless limit 
XSB = NcNf^'^/n^ as /r —oo. There is a clear threshold at a scale /r « Mthr = 0.467GeV. 


Even though it is not strictly proven that a first-order phase transition occurs in cold and dense QCD, it has been 
predicted in several low-energy QGD models. However, no sharp, first-order transition is observed in our predictions. 
The massive parameters in our model, M3 and m are fixed by zero-temperature, zero-^ lattice data for the quark 
propagator. They do not depend on temperature nor chemical potential. A crossover may be reasonably described, but 
a first-order phase transition needs symmetry changes that, we believe, will only be attained by T— and /i—dependent 
parameters or condensates, such as Indeed, M3 can be seen as a nonlocal order parameter for chiral 

symmetry breaking [30j and its T, fi dependence will thus play a significant role in the chiral phase transition. 


C. Results for T A 0 and p 7 ^ 0 


In this subsection, we consider a thermal system with an imbalance between the number of fermions and anti¬ 
fermions, given by a finite /i. In Figs. (6a) and (6b), pressure and quark number density, respectively, are shown as 
functions of chemical potential for different values of temperature. 






Quark Chemical Potential [GeV] 


(a) 


(b) 


FIG. 6: (a) Pressure and (b) quark number density (normalized by the free massless gas result) are plotted as 
functions of the chemical potential fi for temperature T = {0 (bottom), 50,100,150, 200 (top)} MeV. 
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These results further support the physical consistency of our model and the absence of thermodynamical instabilities. 
Moreover, we see that, as temperature increases, more thermal energy becomes available, exciting particles even at 
lower chemical potentials, as expected. Even though we do not see a first-order phase transition at low temperatures, 
as discussed above, we do observe two qualitative features generally seen in QCD models: (i) a smoothening of the 
transition as temperature increases and (ii) a shift of the inflection point of the curve to lower values of /i. 


V. SUMMARY AND DISCUSSION 

In this paper, we made a first exploratory study of the thermodynamics of a softly BRST broken quark model. A 
soft breaking of the BRST symmetry could in principle be present in the infrared regime of QCD without affecting 
any ultraviolet (perturbative) predictions. In this case, the BRST breaking may be intimately related to other 
nonperturbative phenomena, such as confinement and chiral symmetry breaking, representing a complementary way 
to study these issues. Evidence for this breaking has been recently found on lattice studies in the gluon sector [4]. 
It also provides an argument for the absence of asymptotic quarks in the spectrum as well as a good analytical 
description of nonperturbative quark propagators as compared to lattice. 

We adopt the tree-level approximation in our model, so that no explicit quark-gluon interactions are present. 
Nevertheless, gluon dynamics is encoded in the nontrivial BRST-breaking background in the form of a nonlocal quark 
propagator. The model reduces itself to a free quark model in the high energy limit, corresponding to the fermionic 
sector of Quantum Chromodynamics, and yields an infrared dynamics compatible with the available lattice data for 
quark propagators. Within this model, as in other nonperturbative approaches to QCD, confinement is encoded in 
the presence of complex poles in the momentum-space quark propagator, i.e., complex masses. The quark propagator 
adopted here has been explicitly shown to violate reflection positivity, being compatible with the absence of asymptotic 
quark physical states. 

Despite the presence of the nontrivial BRST-breaking background, the theory is quadratic in the quark fields. 
Therefore, the grand canonical partition function of the model could be calculated exactly at finite temperature 
and chemical potential. This allowed a computation of several thermodynamic quantities, as shown in the previous 
section. It is important to note that this model does not lead to thermodynamic instabilities, in spite of the presence 
of complex poles in propagator. We believe that the presence of three poles (one of them real) in the propagator of 
our model is a crucial feature to guarantee the consistency of the thermodynamic predictions, since quark models 
with pair of complex-conjugated poles alone have displayed unphysical macroscopic behavior [42) . 

Besides presenting a consistent thermodynamic behavior - free of instabilities -, the model displays nontrivial 
physical results, even at tree level. At finite temperature and vanishing chemical potential, one sees a nontrivial 
behavior of thermodynamic observables, which clearly differ from the case of free quarks at low and intermediate 
temperatures, while recovering the standard free case at high temperatures. This can be seen as an effect of the 
momentum-dependent mass function, which is in turn a manifestation of the soft BRST breaking. It is interesting to 
notice that there is a clear rise in the pressure, as well as a peak in the trace anomaly, at around T ^ 0.15 GeV, the 
temperature scale of the deconfinement and chiral restoration crossovers. In spite of this, one may not associate this 
change of behavior with a phase transition, as other signatures are not present, as e.g. a dip in the sound velocity. 

The model also allowed us to access the region of finite chemical potential at arbitrary temperatures, including a well- 
defined T —>■ 0 limit. The finite density thermodynamics was also shown to contain nontrivial features, differing from 
an ideal quark gas picture. The system has no excitations until a threshold chemical potential fj, ~ Mfhr = 0.467 GeV 
is achieved, so that the Silver-Blaze property [SH ES] is satisfied. For larger chemical potentials, the limit of massless 
free quarks is approached. In spite of the clear excitation threshold Mthr, it is probably not correct to interpret 
the system as a gas of constituent quarks of mass Mthr, as can be seen, e.g., from the comparison of the respective 
pressures in Fig. It seems more appropriate to resort to an interpretation in terms of an ensemble of resonances 
with different masses. This idea seems to be supported by the momentum and chemical potential dependence of the 
mass function As the chemical potential is increased, more phase space is made available for excitations with 
different momenta. Given that the mass depends on both fj, and p, these newly available states have different masses 
as those present before. A similar interpretation seems to be actually possible for any T or /i. 

Of course, further improvements and perspectives are possible and planned. Our results for the equation of state at 
large densities and low temperatures may in principle be further developed to investigate whether a physical compact 
star with a quark matter core could be created in this model. The explicit inclusion of interactions on top of the 
already nontrivial nonperturbative background may be implemented through different couplings, such as four-fermion 
or quark-gluon vertices, to be added perturbatively. Moreover, in order to explicitly compute the Polyakov loop in 
our model one could compute its effective potential in the background field formalism (for recent pure glue studies, 
see Refs. [551)50] and Ref.[0T| for the case with heavy quarks). This way one may understand further the significant 
change of thermodynamic quantities observed here in terms of the approximate order parameter for the deconfinement 
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phase transition. 
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Appendix A FERMION EFFECTIVE ACTION AT FINITE T AND /r 

A Model hamiltonian 


The first step to calculate the grand partition function is to determine the hamiltonian density of the theory. It 
can be obtained from the lagrangian density in Minkowski space 

U = Coirac + UjA -|- Cm , (49) 

where, omitting spinor, color, and flavor indices, 

Cmrac = 4’ (50) 

and 

C^\ = + m?)X + X{—d^ + — 0{—d^ + m'^)ri, (51) 

and 


Cm = + M2iX4 + ^>^) (52) 

The hamiltonian is the Legendre transformation of the lagrangian with respect to the generalized velocity, that is 

n = Y,^^.i9o4^)-c, (53) 


where II^. is the momentum conjugated to the field (pi- The field momenta are defined as [64] 

dC 


Hci, = 


d{dQpi) 


or 


Ha, = — - 


dC 


for bosonic fields 


for fermionic fields 


(54) 


(55) 


(56) 


It is important to notice that due to the first-order nature of the Dirac equation, the quark fields ip and Tp are actually 
such that iipjo is the momentum conjugated to the field ip-: i-e-j 


didopi) 

After adding a total derivative to C^\, one straightforwardly finds the following set of field momenta 



Bf = -doX 

Ba = doi 

Ba = -do^ 

Be = 

Bg = -^o?? 

B,, = —doB 

II 

tf 


= iip^o 


(57) 


Putting all pieces together, the hamiltonian density can be written as 

^ — ^Dirac + -C 'Em, 


(58) 
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where 


'^Dirac [n^, V’] = -n^7o7 ■ VV' - zmon^7o'!/') 


(59) 


and 


T~L^\ — + n^i9o^ + IIa^o^ + n_;^9oA + Ilg^o^ + + n^9o/y + n^9o7 ~ 

= n^iix + - n^-n^ + n,^ne - (v^) • (va) - (va) • (vo - (v^) • (V0) + (vd) • (Vt?) - 

—m^{^X + X^ + fjO — 0r]), (60) 


and 


'Hm — ~^M — + tp^) + M2{Xip + tpX). 


(61) 


B Noether current 


In the grand canonical ensemble the average charge of the system is kept constant. This constraint can be imposed 
by the replacement "H —> "H — where jo is the charge density and ^ the corresponding chemical potential. From 
a field theoretical point of view, such a charge density is consistently defined as the time component of the Noether 
current associated with a global U{1) transformation. Let us now define one such transformation for the theory at 
hand and find the associated Noether charge. 

In order to define the quark number, we first notice that the action is symmetric under the U{\) transformation 


The corresponding Noether current is 


Ip e-^°‘'ip, 

Ip 





A ^ e-*“A, 

A ^ e“A 


r] ->■ e"“77. 

f] e*“77 


0 e-“6l. 

0 

(62) 

■LL A I 

(63) 


and the charge density is 


jo = X! 


dL 


d{dQ4>i) 






(64) 


/ 


where we recalled the definitions of field momenta for bosons (541 and fermions (551. The explicit expression for the 
charge density in terms of fields and field momenta is 


jo = i'i/'n.,/, + + iAn> — fAn^^ — idllg + i9Ilg — ir]Il^ + irjUfj. 


(65) 


C Quark in-medium effective action 

In the functional integral formalism, the partition function is written as 
Z{T, fj.) =Tr exp 

=J[DU^] [DP;] [D^] [ 17 ^] [DUg] [DX] [DUx] [D~X] [ 170 ;,] [Dr^] [ 770 ,] [Dfj] [DUf,] [D0] [DUe] [D0] [DUg] 

j d'^x \^^mi{dA4>i) -% + fiji 


X exp 


( 66 ) 
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where now the fields have been analytically continued to imaginary time according to xq —t 2^4 = ixQ. The constraint 
of finite temperature is imposed on the system by taking the time direction compact, i.e., 0 < X 4 < j3 = 1/T. This is 
expressed in the integration symbol 



(67) 


Notice that one must start from the expression (661, where integration over field momenta must still be made. As 
is well known from the study of a complex scalar field (see, e.g., [56)1. the field momentum integration in the presence 
of a chemical potential does not simply lead to the lagrangian plus a term. Rather, the result of including a 
chemical potential is typically to shift the time derivative as 84 —>■ (84 — /i). 

The integrand of the exponent of (66) is 


/ ^ iUk^k -l-i + fij4: 

k 

= + illgO + iUgO + til,,?) + + 

+n,/,7o7 ■ V?/' + t?7ion,/,7oV’ - H^IIa - HaII^ + Hgll,, - n,,ne + 

+ (V0 • (VA) + (VA) • (VO + (VO • (V0) - (VO • (V 77 ) + 

+?t?.^(^A + A^ + 7)0 - Orj) + - M 2 {Xtp + ipX) + 

+ i^H^X — 

— + ifiOIlg — + ififjUfj 

= ili^ (^ijj - flip - t7o7 • V?/? + mo-fopj^ + tn^-(l + + *nA(A - ^A) + 

+tnA(A + fiX) + tn^(^ — mO + tng(0 + ^9) + iu^{r] — fj,T]) + 

+?n,,(?7 + fif]) + iiig{9 - fi6) - n^-iiA - + iigii,, - + 

+ (V0 • (VA) + (VA) • (VO + (VO • (V0) - (V0) • (VO + 

+to^(^A + X^ + fjO - Oij) + Ml{£,'ip + ip^) - M2{Xip + ipX) (68) 


where we used the notation {(pk = 84 (pk)- 

Note that we can separetely integrate the bosonic and the fermionic auxiliary fields. As a result, we shall find 
an effective action for the fermion fields. Interestingly enough, we find a similar result as compared to the vacuum 
nonlocal action, Eq. ®, the only difference being the replacement^ 8 \ —)■ —(iujn + /i)^. We now demonstrate this 
result, while we leave the integral in the quark field tp to the main text. 

Let us first consider the anticommuting auxiliary fields (00 A, A). The relevant part of the integrand for the field 
momentum integral is 


Iiif ~II^IIa + *II^(C + /rO + *IIa(A — /iA) — flAlIj + ?nA(A + /rA) + ^0? ~ 


“h % — /^A^ 

+ /iA 

- 

^x + iU- mA 

— i ^A -l- /^A^ 

~ (a — /rA^ — ^A + ^A^ ^ 

1 



Notice that the field variables and momenta above are Grassmann variables. Once this expression is under a functional 
integral, we may shift the momenta in without changing the integral, so that 


Zn, := J[DU^][DU^][DIl4,][Dn-^]efo 

= Jo [(I+m«)(A-?«a)+(a+ma)({-/.{)] y’[ 7 ^n 5 ][i:)n^-][DnA][£>n^]e“-^o'’ d'*x(njnA+nAne) 

_ fo [(f+Aif) (A-MA) + (A+/iA) (5-/i{)j 


^ Notice that this is the same as in the case of Dirac fermions or commuting complex scalar fields in the presence of a chemical potential. 
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Let us now compute the momentum integral of the bosonic auxiliary fields. The relevant integrand is 


iiii = + iUg 


+ iUg ^9 + fj,6^ + zll^ (?7 — /izy) — n^IIe + zll^ {rj + ^rj) + ilVg {o — fiO^ 


= [Ilg + i{r] - ^r])] + z + fj,0^ - - i{9 - ^9) [Ifg - z (zy + /zry)] + 

+ (9 + fi9j (rj - fxr]) - (^ + m) {9 - fx9^ 


(71) 


As in the fermionic case, for fixed field configurations, one is allowed to shift the held momenta by any arbitrary 
function without changing the functional integral. As a result. 


Zn, := J [DIl^][Dn^][Dn0][DUs]e^o 

_ ^lo [(S+^l8'j(r|-^J,r|) + [fj+fJ,fj)(^0-^le)^ 


r* 

[DU^] [DUf,] [Dng] [L'n^Je-^o'’ [nsn„+n,jn<,] 


In our next step, we integrate over the auxiliary boson helds {9,9,r],fi). The relevant integrand is 
h := (j) + ^9^ {fj - fir]) - (?) + fif]) (j) - ]i9^ + (Vzy) • (Vd) - (V0) • (Vry) + m^{f]9 - 9r]) 

= 9f] — fi9ri + ]i9f] — fi^9r] + (V^) • (V6l) + m^f]9 - fjd — ]jLf]9 + ]ifi9 — yz^^zy + (V0) • (Vzy) + zzz^flzy 

—^ 9 {-do + 2]ido - /z^ - + m^) ry - zy {-dg + 2fido - /z^ - + m^) 0 


(72) 


(73) 


where, in the last step we added a total derivative to the integrand. The resulting functional integral is (aside from 
an inhnite numerical constant) 


j[D9][D9][Dfi][Dr]]e^o A = det[-D^ + clet[i^2 _ ^2]-i(4ZV.) 

= detl-D"^+ (74) 

where we dehned the operator := (^4 — yz)^ + V^, which corresponds to (zw„ + yz)^ — in momentum space. At 
this point, it should be clear that the hnal result should be the same as the one at zero temperature (in euclidean 
space), with the replacement p4 —>■ a;„ — zyz. Let us continue our calculation in order to explicitly check that this is 
the case. The relevant integrand in the fermionic auxiliary helds is 

If := (I + fit) (a - m) + (U fiX) (i - fit) + (VO • (VA) + (VA) • (VO + mO^A + AO + 

+ z/^O - M 2 {X'tfj + rpX) 

^ 5 {-D"^ + m2) A + A (-7?2 + m2) e + z/^O - + z/^A), (75) 


where we added a total derivative and used the dehnition of in the last step. 

The resulting functional integral is again not altered if we add constant functions to the integrated helds. Once the 
z/) held is kept hxed in the integration of (00 A, A), we may freely make the shifts 

^ ^ C + A^2A-V C ^ ? + M2Z/;A-1 
A ^ A-MfA-V A ^ A-M20 A-i 

(76) 


where A := —D^ + m?. It follows that 


0 = (^ + M2 z/>A-i) a (a - m2A-V) + (A - M2z/;A-i) A (^ + MaA"V) + 

+M2 {^ + M2ipA~^) ^|J + Mfip + M2 A~^z/)) — M 2 (a — M'^rpA~^) ^|J — M2'ip (A — M^A^^z/;) 

= CAA + AA^ + 2M2M2z/;A-V (77) 


The integral in the fermionic auxiliary helds is then 


J [D^][D^][DX][DX]e- ^0 


det[A]^^‘=e"'^o'' 

det[-7l2 +^2]4ZV.g-/'’ d-x2MlM2^{^^)i, 


(78) 
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We now finally arrive at our expression for the partition function in terms of the quark fields only 

rP 


Z{T,^i) = J [DIl^][D 2 p]exp 
where the nonlocal mass operator is 


J d^x 111^74 ( 74(^4 - Ai) - Z7 • V + Mj^ Ip 


Mr,= 


2MfM2 


Too 


2MIM2 


mo 


Finally, after the change of variables'^ 11^ —>■ = — tn,^74, the partition function can be written as 




where 


= i’ 


d^xCnl [ipA] 


2MfM2 


74(94 - -\ - \ \ + mo 


iP 


(79) 


(80) 


(81) 


(82) 


is the (nonlocal) quark effective lagrangian at finite temperature and chemical potential. 
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